Limits on cosmological variation of quark masses and strong interaction 
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We discuss limits on variation of {mq/ Aqcd)- The results are obtained by studying n — a- 
interaction during Big Bang, Oklo natural nuclear reactor data and limits on varition of the proton 
g-factor from quasar absorpion spectra. 
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I. INTRODUCTION 

Recent astronomical data suggest a possible variation 
of the fine structure constant a = /Tic at the 10~^ level 
over a time-scale of 10 billion years, see Q (a discussion 
of other limits can be found in Ref. ||] and references 
therein). The data motivated immediately more general 
discussions of possible variations of other constants. Un- 
like for the electroweak forces for the strong interaction 
there is generally no direct relation between the coupling 
constants and observable quantities. In a recent paper 
we presented general discussions of possible influence 
of the strong scale variation on primordial Big Bang Nu- 
cleosynthesis (BBN) yields, the Oklo natural nuclear re- 
actor, quasar absorption spectra and atomic clocks. 

Here we continue this work, concentrating on the prop- 
erties of ^He (BBN), Sm (Oklo natural nuclear reactor), 
^^C (stars) and proton magnetic g- factor (quasar absorp- 
tion spectra). Since one can measure only variation of the 
dimentionless quantities, we can extract from our results 
variation of the dimensionless ratio mq/ Kqcd where mq 
is the quark mass (with the dependence on the normaliza- 
tion point removed) and Aqcd is the strong interaction 
scale determined by a position of the pole in perturba- 
tive QCD running coupling constant. It is convinient to 
assume that Aqcd is constant, and quark mass mq is 
variable. 

A position of an energy level in a strong potential de- 
pends both on the parameters of the strong interaction 
and nucleon mass M. For example, in a simplest case of a 
particle in a deep square potential well the single-particle 
energy levels are 
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(1) 



where V and R are the depth and radius of the po- 
tential. The strong interaction is sensitive to the light 
quark mass because the 7r-meson mass m,r ~ ^^mqAQCD, 



where m„ = 



-md. 



The nucleon mass is also sensitive 



to quark mass because M has large "strange fraction" Q| 



dM 
dmg 



Putting the strange quark mass TOs = 130 MeV one 
finds that about 1/5 of the nucleon mass comes from 
the "strange term" . 



II. N - a INTERACTION 

The production of nuclei with A > b during BBN is 
strongly suppressed because of the absence of stable nu- 
clei with A = h. ^He is unstable nucleus which is seen 
as a resonance in n-a elastic scattering at neutron labo- 
ratory energy around 1.1 MeV The resonance corre- 
sponds to the ground state of unstable nucleus '^He with 
mass excess A = 11386.234 KeV The ground state 
lies at 0.89 MeV above neutron threshold. This energy 
is rather small in nuclear scale, therefore, rather small 
variation of the strong interaction scale may influence 
considerably the position of the resonance making, for ex- 
ample, ^He stable. Stable ^He at the time of BBN would 
change strongly the primordial abundances of light ele- 
ments. Basing on the fact that the standard BBN theory 
describes rather well the observed light element abun- 
dances we can put bounds on strong scale variation at 
the time of BBN. 

A phenomenological potential describing n — a scat- 
tering phases at low energy has the Woods-Saxon shape 

i- 

U{r) = -T/(e^ + l) + (n/™^c)V,(L-cr)r-i(d/dr)(e^=-|-l). 

(3) 



In Eq.i), 



x = [r-R)/a, Xs^{r-Rs)/as 



(4) 



=< N\ss\N >« 1.5 



(2) 



Here, R and Rs are the central and the spin-orbit radiuses 
of the potential, 

R = roiMjMpy/\ Rs - r,(M,/Afp)i/3, 

where (Ma/Mp) = 3.973 is the ratio of the a and the 
proton masses. The other parameters of the potential 
are 

y = 41.88 MeV, K = (3.0 + 0.1£:„) MeV 

ro = (1.50 - 0.01£:„) fm, = 1.0 fm 
a = fls = 0.25 fm 

The main contributions to this potential come from ex- 
change of scalar and vector mesons. These contributions 
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are not related directly to the pion and we can expect 
them to be insensitive to variation of the pion proper- 
ties. Pion exchange interaction gives zero contribution 
into the potential via direct term. The exchange term, 
however, gives nonzero contribution. If we isolate this 
contribution, we can study, at least approximately, the 
dependence of ^He ground state properties on pion mass 
and pion-nucleon coupling constant. 



A. One pion exchange (OPE) contribution 



position E = 0.89 MeV. The procedure gives V = 32.74 
MeV that should be compared with 41.88 MeV for full 
phenomenological potential Eq.(||). The difference be- 
tween these two values can be treated as an effective 
depth of equivalent local potential corresponding to non- 
local exchange contribution of pion, Eq.(|). The Eq.(|) 
was obtained neglecting pion-nucleon formfactor. This is 
reasonable approximation since the cutoff parameter A in 
the formfactor is large compared to nucleon momentum 
in a-particle nucleus. 



We start from the standard OPE interaction 
P 



K(ri-r2) 



2 {Tl-T2){fT, ■ Vi)(<T2 • V2)e-'"''--/ri2. 



(5) 

Here we distinguish between pion mass in the exponent 
and pion mass in the denominator. The last one is used 
for normalization of pion-nucleon coupling constant and 
should not be varied. 

The exchange contribution into n — a potential can be 
written in the form 



47r fi^ ri2 



^i?o(r-2). (6) 



In Eq.(||), Ro{r) is Is radial wave function of nucleon 
inside a-particle. For our purpose it is sufficient to choose 
Ra{r) as a simple gaussian 

Roir)^^ily/''expi-ry'ib^), (7) 

where the parameter b is related to mean-square size of 
the a-particle 

We omitted the contact term in Eq.(^. It does not de- 
pend on the pion mass and can be included into phe- 
nomenological part of the potential Eq. (||) together with 
other mesons contribution. With the potential Eq.(||) 
the Schrodinger equation for n — a system becomes non- 
local. For P3/2 radial wave function we obtain the 
equation 
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Xi{r')dr' 
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In Eq.(g), ii{x), and ki{x) are the spherical Bessel func- 
tions of imaginary argument. 

Using Eq.(||) we can now adjust the phenomenological 
parameter V putting the resonance in P3/2 wave in its 



B. BBN bounds on pion, nucleon and quark mass 
variation 

Using Eq. (||) wc can study the dependence of /2 res- 
onance position on pion mass. The potential Eq.(|^) is 
proportional to pion mass squared. Therefore, an in- 
crease in pion mass would lead to lowering of the reso- 
nance position. Simple numerical exercise leads to con- 
clusion that the increase in pion mass on 21.9% would 
produce ^He nucleus with zero neutron binding energy. 
Any further increase of pion mass would lead to stable 
^He producing serious modifications in the process of nu- 
cleosynthesis. Using the relation for the vr-meson mass 



\Jvriqh.QCD we obtain the limit 



5{mq/KQCD) 
(mq/AQCD) 



< 0.4 



(9) 



Note that this OPE path to obtain the limits does not 
look very reliable - see discussions in recent papers 
0. Therefore, below we present the limits obtained in a 
different way. 

Variation of the nucleon mass M leads to change of 
the kinetic energy and change of the resonance position. 
We found that 5% increase of the nucleon mass produces 
^He bound state with zero neutron binding energy. Using 
the relation (H) we find the limit for the variation of the 
strange quark mass 



6{ms/AQCD) 
{nis/AQCD) 



< 0.25 



(10) 



These are limits on the variations of the light quark mass 
and strange quark mass from BBN to present time. 



III. LIMITS FROM OKLO NATURAL 
NUCLEAR REACTOR 

Similar limits follow from data on natural nuclear re- 
actor in Oklo active about 2 bn years ago. The most 
sensitive phenomenon (used previously for limits on the 
variations of the electromagnetic a) is disappearance of 
certain isotopes (especially ^^^Sm) possessing a neutron 
resonance close to zero Today the lowest resonance 
energy Eq — 0.0973 ± 0.0002 eV is larger compared to its 
width, so the neutron capture cross section a ~ ^/Eq- 
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The data constrain the ratio of this cross section to the 
non-resonance one which was used to find neutron flux. It 
therefore imphes ( under assumption that the same reso- 
nance was the lowest one at the time of Oklo reactor) that 
these data constrain the variation of the following ratio 
5{Eq/Ei) where Ei ^ IMeV is a typical single-particle 
energy scale which is determined mainly by Aqcd- 

A very small value of Eq appears as a result of the 
nearly exact cancellation of two large terms Eq — Ec — 5„ 
where Sn is the neutron separation energy in ^^''Sm and 
Ec is the energy of the many-body excited compound 
state relative to the ground state of ^^°Sm (recall that we 



^He calculation (accuracy in this term is very low any- 
way). 

The limit on the shift of the Sm resonance is ^ \SEo\ < 
0.02 eV. A comparison with Eq. (|l5|) gives very stringent 
limits: 



consider the reaction 
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Sm) . The difference is very 



small because we deliberately selected the lowest com- 
pound resonance above neutron threshold. The neutron 
separation energy is Sn = V — ep where V is the depth 
of the potential and ep = p)/2M - {Ti^ A^/^)/{MR^) is 
the Fermi energy. The energy of compound state Ec is 
approximately equal to the difference of several single- 
particle energies given by Eq. (|^). Therefore, both ep 
and Ec scale as Tt?/{MB?). We can present the reso- 
nance energy in the forllowing form: 



Eo=Ec-Sn = Ec + ep-V = K 



V (11) 



where _ftr is a numerical constant. We can find this con- 
stant from the present time condition Eq w 0: 



K {, pre sent — ^present 



(12) 



Now we may consider change of Eq under variation of 
M, i?, V. 



Comparing with the ^He calculation in the previous sec- 
tion we can say that two last terms in the brackets cor- 
respond to the variation of the energy level due to the 
variation of the pion mass m^. This gives us the follow- 
ing expression 

dE ^ dE dE ^ 

SEq = -V— + ^Sm. = —SM + -^Sm. (14) 

We can find numerical values of the derivatives using V— 
50 MeV, SM/M = 0.2(5m^/m^ (see Eq. (|)), Sm^/m^ = 
0.5dmq/mq. 



SEq 



-0.05JM-0.03(5m^ = -10 MeV 



SrUs 



-2MeV^-^ 
ruq 

(15) 

Note that in the ^He calculation the derivative 4^ had 
smaller absolute value (-0.02 instead of -0.05 here). The 
difference may be related to the fact that the p-wave reso- 
nance wave function in ^He was localised mainly outside 
the nucleus (where the coefficient V in Eq. (p^) is equal 
to zero). We extrapolated the derivative ^^^^from the 



.S{M/Aqcd) , ^Sim^/AQco) 



(M/Aqcd) 



0.6- 



(m^/Agcn) 



< 4 X 10" 
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S{ms/AQCD) , „ „<^(™<;/Aqcd) 



(TOs/Aqcd) 



0.2- 



(mq/AQCD) 



< 2 X 10" 



(16) 



(17) 



Note that the authors of the last work in [gj found also 
the non-zero solution SEq = -0.097 ± 0.008 eV. This 
solution corresponds to the same resonance moved below 
thermal neutron energy. In this case 



S{M/Aqcd) , ^^Sim^/AQCD) 



(M/Aqcd) 



(m^/AgcD) 



= (2±0.2)xl0"^ (18) 



S{ms/AQCD) , Q ^5{mq/AQCD) 



{ms/AQCD) 



{mq/AQco) 



= (l±0.1)xlO-^ (19) 



In principle, the total number of the solutions can be very 
large since ^^"Sm nucleus has millions of resonances and 
each of them can provide two new solutions (thermal neu- 
tron energy on the right tale or left tale of the resonance). 
However, these extra solutions are probably excluded by 
the measurements of the neutron capture cross-sections 
for other nuclei since no significant changes have been 
observed there also, see H. 



r ,SM 2SR, , SM , 2SR , , 

= -^MR^^W^^^-'"" = 

(13) 



IV. LIMITS FROM ^^C PRODUCTION IN 
STARS AND QUASAR ABSORPTION SPECTRA 

In the previous sections we obtained limits on variation 
of rris/ Aqcd during the interval between the Big Bang 
and present time and on shorter time scale from Oklo 
natural nuclear reactor which was active 1.8 billion years 
ago. It is also possible to obtain limits on the interme- 
diate time scale. One possibility is related to postion of 
the resonance in ^^C during production of this element 
in stars. This famous resonance at E=380 KeV is needed 
to produce enough carbon and create life. According to 
Ref. 1^ the position of this resonance can not shift by 
more than 60 KeV ( one can also find in Ref. |^ the 
limits on the strong interactions and other relevant ref- 
erences). We can use Eq. (^sj) to provide some rough 
estimates only: 



S{M/Aqcd) , ^J{m^/AQCD) 



{M/Aqcd) 



0.6- 



{m^/AQCD) 



< 1.2 X 10" 



<?(m /Agcz.) ^Q /(m /A ) ^ ^ ^ ^^^3 
(ms/AQCDj [mq/AQCD) 



(20) 



(21) 



4 



Unfortunetely, these limits are relatively weak. Much 
stronger limits can be obtained from the measurements 
of quasar absorption spectra. Comparison of atomic H 
21 cm (hyperfine) transition with molecular rotational 
transitions gave the following limits on Y = a^gp | |lO[ | 
SY/Y = (-.20 ± 0.44)10-5 for redshift z=0.2467 and 
SY/Y = (-.16 ± 0.54)10-5 for z=0.6847. The second 
limit corresponds to roughly t=6 billion years ago. These 
results were used in to obtain the limits on variation 
of a. In Ref. Q it was suggested to use SY/Y to estimate 
variation of mg/ Aqcd- Here we continue this work. 

According to calculation in Ref. |ll| dependence of 
the proton magnetic moment on 7r-meson mass can 
be approximated by the following equation 



Mp(0) 



1 



(22) 



where a= 1.37/GeV,h^ 0.452/GeV2. The proton g- 
factor is gp = fj,p/{eh/2Mc). Using Eq. (H) we obtain 
the following estimate: 

Sor, SM SttItt (5™, 6m„ 

^ 0.18 ^=0.2 ^-0.09 ^ (23) 

gp M TOtt ms ruq 

Comparing with the limits on SY/Y we obtain for the 
red shift z=0.2467 (neglecting variation of a): 

5{M/Kqcd) _Q^^^5{mJKQCD) ^ (_o.2±0.44) x 10"^ 



{M/Kqcd) ■ {m^/^QCo) 

^ SjiTis/AQCD) _ S{mq/AQCD) 
(ms/Agcu) (mg/Agcn) 



(24) 



(-2.0 ±4.4) X 10" 



For the red shift z=0.6847 : 

S{M/Aqcd) q ^ Jim^/AQCP) 
(M/Aqcd) ■ (tox/Aqcd) 

^ S{ms/AQCD) _ 6{mq/AQCD) 
{itIs/Aqcd) {mq/AqcD) 



= (-0.16±0.54)xl0-5 
(26) 

= (-1.6 ±5.4) x 10-5 



_ (27) 

There is also a limit on A = a^gpirie/Mp |1^ SX/X = 
(0.7 ± 1.1)10-5 for z=1.8. This limit was interpreted as 
a limit on variation of a or rUe/Mp. 

Using eq. (p2|) we obtain 



SigpmjM) _ Sm^ ^ 18^^ ^ ^ _o 09^ (28) 



me 



(gpUie/M) m. 
This gives us a limit for the red shift z=1.8 



SimJAQcn) ^ ^/(m./AQcz.) ^ ^^.y ± 1.1) x 10-5 



(TOe/Agci)) 



{mq/AQco) 



(25) 



(29) 

We should stress that in this paper we present experi- 
mental errors only. The theoretical errors are hard to 
estimate, therefore some of these limits, possibly, should 
be treated as order of magnitude estimates. 
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